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ABSTRACT 


To  Investigate  the  possibility  of  a point  of  a crack  front 
propagating  In  a direction  not  In  the  normal  plane  of  the  crack 
front,  the  three-dimensional  form  of  the  crack  front  stress 
field  Is  obtained.  To  simplify  comparisons  of  the  states  of 
stress  at  various  points  lying  on  a spherical  surface  centered 
at  a point  of  the  crack  front,  a local  spherical  coordinate 
system  Is  used.  It  Is  found  that  the  crack  propagation  will 
be  from  each  point  of  the  crack  front  In  a direction  lying  In 
the  normal  plane. 

The  results  are  used  In  conjunction  with  the  strain  energy 
density  fracture  criterion  for  the  problem  of  an  elliptical 
crack.  The  plane  of  the  flat  elliptical  crack  makes  an  arbi- 
trary angle  with  the  field  of  uniform  applied  tensile  stress. 
Crack  growth  directions  for  various  positions  along  the  crack 
front  are  determined,  and  loads  required  for  fracture  for  var- 
ious angles  are  obtained. 


INTRODUCTION 


Cracks  of  particular  shapes  have  been  considered  by  sev> 
eral  authors.  The  plane  strain  geometry  and  loading  always 
result  In  the  familiar  plane  strain  crack  tip  stress  singular- 
Ity  established  by  Williams  [1].  Similarly,  the  plane  strain 
geometry  subjected  to  anti-plane  loading  always  leads  to  the 
form  of  crack  tip  stress  singularity  given  In  [2].  These  forms 
win  be  referred  to  as  the  two-dimensional  results.  They  con- 
tain the  three  stress  Intensity  factors,  k^ , V.2*  and  k^.  Ex- 
amination of  the  two-dimensional  results  shows  that  the  singu- 
lar part  of  the  crack  tip  stresses  depends  explicitly  on  posi- 
tion relative  to  the  crack,  and  Implicitly  on  the  other  param- 
eters (load,  geometry,  material)  through  the  stress  intensity 
factors. 

This  same  type  of  two-dimensional  geometry  subjected  to 
arbitrary  three-dimensional  loading  was  the  subject  of  an 
eigenfunction  expansion  Investigation  [3].  It  was  found  that 
the  crack  tip  stress  singularity  was  Identical  In  form  to  the 
two-dimensional  results,  with  the  exception  that  the  stress 
Intensity  factors  depend.  In  general,  on  the  position  of  the 
point  of  the  crack  front  nearest  the  point  at  which  the  stresses 
are  computed. 

Irwin  £4]  showed  that  the  form  taken  by  the  stress  field 
in  [5]  for  a penny  shaped  crack  Is  also  the  same  as  the  two- 
dimensional  result.  Progress  since  then  has  established  that 
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this  Is  true  for  very  general  loads.  For  a review  of  the 
problem  of  the  penny  shaped  crack,  see  [6],  which  contains  a 
large  number  of  stress  Intensity  factors  evaluated  for  vari- 
ous loads. 

Finally,  for  the  cases  considered,  the  form  of  the  stress 
field  near  the  edge  of  an  elliptical  crack  is  also  given  by 
the  two-dimensional  results  £7].  Reference  [6]  contains  sev- 
eral examples  of  loading  and  the  corresponding  stress  inten- 
sity factors. 

Although  not  proved,  it  seems  quite  likely  that  similar 
results  apply  to  arbitrary  cracks.  The  starting  point  in 
this  investigation  is  the  assumption  that  the  two-dimensional 
results  give  the  form  of  the  stress  singularity  for  an  arbi- 
trary crack.  The  crack  front  is  a general  curve  in  space, 
described  implicitly  by  its  curvature  and  torsion.  It  will 
also  be  assumed  that  the  osculating  plane  of  the  crack  front 
curve  is  tangent  to  the  free  surfaces  of  the  crack  at  each 
point. 

A geometrical  analysis  provides  the  stress  field  referred 
to  local  spherical  coordinates.  The  strain  energy  density 
theory  [8]  when  applied  In  three-dimensions  [9]  requires  this 
result  so  that  minimums  of  strain  energy  density  on  spherical 
surfaces  can  be  located.  The  result  is  surprisingly  simple 
when  the  appropriate  choice  of  coordinate  system  is  made.  A 
different  choice  of  coordinates  can  lead  to  results  as  compli- 
cated as  those  in  Appendix  I £6,9]. 
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The  results  are  applied  to  mixed  mode  loading  of  a flat 
elliptical  crack.  The  stress  Intensity  factors  of  [7]  can 
be  combined  to  give  the  crack  front  stress  singularity  for  a 
field  of  uniform  tension  at  an  arbitrary  angle  to  the  plane 
of  the  crack.  The  various  crack  growth  directions  and  loads 
are  displayed  In  Appendix  II  as  functions  of  the  parameters 
1 nvol ved . 

A future  application  of  this  result  to  the  transverse 
crack  In  a thick  plate  Is  anticipated.  The  development  of  a 
shear  lip  due  to  the  Interaction  of  a curved  (thumbnail)  crack 
front  with  the  free  surface  will  lead  to  mixed  mode  crack  front 
stress  fields.  The  analysis  of  additional  crack  growth  under 
mixed  mode  loading  requires  further  work  on  the  finite  element 
numerical  stress  analysis.  But  when  this  Is  obtained,  the 
strain  energy  density  fracture  criterion  will  be  available  for 
direct  application. 

Further  work  to  Include  higher  order  terms  seems  appro- 
priate for  use  In  developing  special  crack  front  elements  In 
finite  element  formulations  of  problems  Involving  curved 
cracks.  It  Is  expected  that  the  curvature  and  torsion  of  the 
curve  and  their  derivatives  along  the  crack  will  be  required. 
Derivatives  of  the  stress  Intensity  factors  may  also  appear 
In  the  result. 


-4- 


LOCAL  STRESS  FIELD 


Let  P be  a point  of  the  curve  defining  the  crack  front  and 
let  n,  b,  and  t be  the  normal,  blnormal,  and  tangent  Cunit) 
vectors  at  P.  The  plane  defined  by  the  normal  and  tangent, 
known  as  the  osculating  plane,  will  be  assumed  to  coincide 
with  the  tangent  to  the  free  surfaces  of  the  crack  at  P.  This 
Is  by  no  means  necessary,  but  for  a naturally  growing  crack. 

It  seems  reasonable.  For  an  element  In  the  normal  plane,  de- 
fined by  n^  and  b.  as  shown  In  Figure  1.  the  appropriate  form 
of  the  stress  field  Is  the  two-dimensional  result. 
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Figure  1:  Crack  front  stresses  In  the  normal  plane. 


These  stresses  act  on  an  element  with  edges  In  the  n,  b.  and 
^ directions.  The  stresses  on  the  p1ane»  9»0  reduce  to 


In  all  cases,  , k2>  and  k^  are  to  be  evaluated  at  point  P 
of  the  crack  front.  The  expressions  for  the  stresses  oro1t 
terms  of  0(1).  which  are  finite  at  the  crack  front. 


GEOMETRY 


Consider  two  points,  and  P,  on  the  crack  front  as 
shown  In  Figure  2.  Regard  P^  as  fixed  and  P as  moving.  Then 


Figure  2:  Vector  triads  at  two  points. 

the  position  of  P,  rp,  and  the  unit  vectors,  n,  b,  and  t may 
be  regarded  as  functions  of  the  arc  length,  s,  from  P^  to  P. 
It  Is  shown  In  differential  geometry  [10]  that 
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dt 

H?  “ ‘ "5 

H?  " - *'0 

The  sign  of  the  curvature,  ic.  Is  opposite  that  In  [10]  because 
n here  points  to  the  convex  side  of  the  curve  when  k Is  posi- 
tive. When  the  torsion,  t.  Is  positive,  n and  b turn  In  the 
direction  of  a right-handed  screw  about  t as  P moves  In  the 
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direction  of  t, 


If  s is  small,  Taylor  series  expansions  of  the  vectors 
about  S“0  are 
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where  equations  (3)  have  been  used.  The  zero  subscript  in- 
dicates that  the  term  is  evaluated  at  P . 

0 

The  position  of  an  arbitrary  point  Q may  be  expressed  in 

two  ways.  In  terms  of  its  coordinates,  n^^,  b^,  and  t^^,  in 

the  reference  frame  at  P , 

0 


To  “ »'p  + b„b„  + t„t^ 

" ^0  ®"® 


Using  the  coordinates,  n,  b,  and  t,  in  the  reference  frame 
at  P, 


Tq  ■ rp  + nn  + bb  + tt 


r 

I 


Equating  (5)  and  (6)  and  using  equations  (4)  the  relation  be> 
tween  the  two  sets  of  coordinates  Is  found  to  be 


”o  " " “ ^*^0^  * 

b©  = b + Tj^ns  + 0(s*)  (7) 

to  = t + (1  + 

The  direction  cosines  needed  for  transforming  the  stress 
components  from  one  reference  frame  to  the  other  may  be  ex- 
pressed as  In  Table  I.  The  elements  In  the  table  are  cosines 


Table  I:  Direction  cosines  of  the  transformation 


of  the  angles  between  the  two  sets  of  unit  vectors.  These 
are  easily  obtained  from  equations  (4).  Note  that  terms  of 
O(s^)  have  been  omitted  from  Table  I. 
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TRANSFORMATION  OF  COORDINATES 


Suppose  the  point  Q of  Figure  2 is  the  location  of  the 
elements  shown  in  Figures  1 and  3.  From  Figures  1 and  3,  the 
coordinates  of  Q,  which  are  related  by  equations  (7),  may  be 


Figure  3:  Crack  front  stresses  in  an  arbitrary 

position. 

written  as 

n * r COS0  n„  = r„  cosip^  cos6„ 

0 0 0 0 

b = r sine  b„  = r„  cos^^  sin0„  (8) 

t - 0 ^0  “ ''o  ^^"*0 

Then  equations  (7)  become 

-10- 


Tq  coS(J>g  COS0Q  = r cose  - sine  s + O(s^) 

Tq  cos4>q  sine^  = r sine  + x^r  cose  s + O(s^)  (9) 

To  sin4)^  = (1  + cose)s  + O(s^) 

For  the  present  investigation,  and  the  spherical  co- 
ordinates, r^,  e^,  and  4)^,  are  regarded  as  given.  Equations 
(9)  are  then  used  to  determine  r,  e,  and  s.  That  is,  the  po- 
sition of  the  point  P such  that  Q lies  in  the  normal  plane  at 
P is  to  be  determined.  The  first  two  of  equations  (9)  give 

r sine  = r^  cos^i^  [sine^  - x^^  cose^^  s]  + 0(s*) 

(10) 

r cose  * r^j  cos4|q  [cose^  + x^  sine^^  s]  + O(s^) 
and  the  third  becomes 

<"o  sin4>jj  = (1  + K^r^j  cos4>q  cos6q)s  + 0(sM  (11) 

Equation  (11)  has  the  solution 

s = Tq  sin<|>Q  + 0(rj)  (12) 

and  equations  (10)  become 


« 


(13) 


r sine  = cos(|i„  sin0^  + O(r^) 

0 0 0 0 

r cose  = r„  cos4)„  cose^  + 0(rM 
0 0 0 0 

In  view  of  equations  (13),  the  coordinates  r and  e are 
r = r^  coS(J)jj  + 0(rj) 


e = 6„  + 0(r„) 


From  these,  one  obtains  such  results  as 


/27  /2r„cos(|i^ 

0 0 


[1  + 0(rj] 


sine  = slne^  + 0(rg),  cose  » cose^^  + O(r^) 


(14) 


(15) 


for  use  In  equations  (1). 


TRANSFORMATION  OF  STRESS 

The  two  sets  of  stress  components  shown  In  Figures  1 and 
3 both  describe  the  state  of  stress  at  the  same  point  Q. 

Thus  they  are  related  by  the  standard  transformation  formulas 
for  stress.  For  example. 
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n nn  n no  b nt  t nn  nb  nb 


^ ^^nn^nt'^nt  * ^^nb^^nt'^bt 


■^nb  “ ®'nn*'bn'^n  * ^nb^bb°b  * ^'nt^bt'^t 


(16) 


■*■  ^^'nn^bb'^nb  ■*■  ^*’nn^bt'^nt  * ^*'nb^bt'^bt 


where  £,u  = a^‘b  is  the  cosine  of  the  angle  between  the  a„ 
direction  and  the  b direction.  These  are  shown  in  Table  I, 
where  the  first  subscript  on  indicates  the  proper  column 
and  the  second,  the  row. 

The  result  for  a°  is  given  to  show  the  order  of  approxi- 
mation involved  in  the  stress  transformation  due  to  neglecting 
terms  in  the  direction  cosines. 


= o 
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- 


* Vnb 


)s  + 0(s* ) 


(17) 


Similar  results  can  be  written  for  the  other  stress  compo- 
nents, but  since  s is  given  by  equation  (12),  the  linear  terms 
In  s are  of  the  same  order  of  magnitude  as  the  terms  neglected 
in  equations  (1).  Therefore,  to  the  order  of  accuracy  required. 


oj  - On  + 0(s) 
aj  » Ojj  + 0(s) 

oj  ■ Of 


- Tnb  + 0(s) 
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LOCAL  STRESS  FIELD  IN  SPHERICAL  COORDINATES 


The  stress  components  in  Figure  3 can  be  evaluated  after 
one  additional  consideration.  The  stress  intensity  factors 
in  equations  (1)  are  to  be  evaluated  at  point  P.  However. 
Taylor  series  expansions  can  be  used  to  relate  these  values 
to  the  values  at  point  P^.  Thus, 

= k°  + 0(s) 

4 * k2  + 0(s)  (19) 

k3  = k^  + O(s) 


By  assembling  together  equations  (1),  (12),  (14),  (15), 
(18),  and  (19),  the  desired  result  is  obtained. 
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ol  = 2v  [— 

/2r.cos4> 


cos  ^ 


/2r„cos^„ 
0 0 


sin  ^3 


nb 


/2r^cos(>^ 
0 0 


0 e 30 

s1n  2~  cos  ^ cos 


/2r~coT^ 
0 0 


cos 


6«  30„ 

2~  E ^ “ sin  -j—  sin  — ^ 3 


(20) 


^ Sin  / 


/2r^cos<))^ 
0 0 


0 •'3  6o 

“ - cos  ^ 

^ /2rcos<t>  ^ 

0 0 


In  equations  (20),  terms  which  remain  finite  as  r^j-^O  are 
omitted. 


Thus,  the  spherical  coordinates  of  Figure  3 result  in  a 
simple  form  for  the  crack  front  stress  singularity.  It  should 
be  noted  that,  to  the  order  of  approximation  used,  the  stresses 
on  the  element  of  Figure  3 are  the  same  as  those  on  an  element 
in  the  normal  (Dq^o^  plane  located  a distance  r^  cos^p  from 
Pjj  and  an  angle  0^  from  n^^.  The  more  complicated  form  re- 
sulting from  the  choice  of  coordinates  used  in  [6,93  is  given 
in  Appendix  I. 

STRAIN  ENERGY  DENSITY 

For  plane  problems,  the  strain  energy  density  fracture 

criterion  [83  states  that  crack  growth  is  toward  the  point 
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I 

I 

on  a circle  centered  at  the  crack  tip  which  has  the  smallest 
I strain  energy  density.  Furthermore,  the  growth  begins  when 

this  smallest  value  of  the  strain  energy  density  reaches  the 
i allowable  value  for  the  material. 

( A similar  criterion  has  been  stated  [9]  for  three-dimen- 

sional problems.  Each  point  of  the  crack  front  must  be  con- 
* sidered  separately.  A small  sphere,  centered  at  a particular 

j point  of  the  crack  front,  is  used  to  replace  the  circle  In 

the  two-dimensional  case.  The  criterion  states  that  crack 
I growth  is  directed  along  the  line  from  the  center  of  the 

sphere  to  the  point  on  the  sphere  with  the  smallest  strain 
‘ energy  density.  The  collection  of  these  lines,  one  for  each 

I point  of  the  crack  front,  defines  the  new  crack  surface. 

Growth  along  these  directions  begins  when  the  value  of  mini- 
( mum  strain  energy  density  reaches  the  maximum  which  the  ma- 

terial will  tolerate. 

In  use,  since  some  points  fail  before  others,  it  is  con- 
venient to  imagine  an  increment  of  crack  growth  proportional 
in  length  to  the  value  of  minimum  strain  energy  density  at 
each  point.  Thus,  the  points  failing  first  have  more  new 
crack  surface  adjacent  to  them.  The  usual  warning  £ll,12j 
that  the  calculations  apply  only  to  the  first  crack  growth 
from  the  initial  geometry  is  appropriate.  For  very  brittle 
^ materials,  it  has  been  noted  [11]  that  predictions  based  on 

the  initial  crack  geometry  agree  well  with  experiment  for 

large  amounts  of  crack  growth. 
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Since  the  form  of  equations  (20)  Is  so  much  like  the 
two-dimensional  result,  much  use  can  be  made  of  I8J.  The 
strain  energy  density  can,  in  fact,  be  written  easily  as 


^ - S 

dV  " r„cos^^ 
0 0 

where 


+ 0(1) 


(21) 


S = a^^Ck?)"  + 2a^2*'l''2  ^ + 833^^5)"  (22) 

The  coefficients  in  equation  (22)  are  exactly  the  same  as  in 

[8], 


16  Ga^^  * (3  - 4v  - cose^)(l  + cosO^) 

16  0a^2  * 2sin0Q  - 1 + 2v) 

(23) 

16  Ga22  * 4(l-v)(l  - cose^)  + (3cos0p  - 1)(1  + cos0^) 

16  Ga33  * ^ 

where  6 and  v are  the  shear  modulus  and  Poisson’s  ratio,  re- 
spectively. 

When  dependence  on  the  angle  Is  analyzed,  it  is  seen 
that  the  minimum  of  equation  (21)  occurs  in  the  normal  plane, 

- 0 (24) 
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I 

I 
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The  remainder  of  the  analysis  Is  similar  In  some  ways  to 

that  In  [8].  Given  the  stress  Intensity  factors,  one  computes 

the  angle,  9*  which  minimizes  equation  (22).  The  minimum 

* 

value  Is  denoted  by  S * Then  crack  growth  occurs  when 

S*  ■ (25) 

4r  it 

In  the  direction  defined  by  = 0 and  0^  * 6^.  The 

0 0 0 0 

plane  strain  fracture  toughness  denoted  here  by  Is  equal 
to  that  defined  by  ASTM  divided  by  /n. 

The  results  differ  from  two-dimensional  analyses  by  vir- 
tue of  the  fact  that  the  stress  Intensity  factors,  In  general, 
vary  from  point  to  point  along  the  crack  front.  Thus,  the 
new  crack  surface  Is  generally  a warped  surface  of  the  type 
shown  in  (6,9]. 

The  problem  of  uniform  tension  of  an  embedded  elliptical 
crack  Is  presented  as  an  example  In  Appendix  II.  This  prob- 
lem has  been  discussed  In  [6,93,  but  Appendix  II  contains 
several  results  correcting  the  details  of  [6,9]. 
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APPENDIX  I 


I 

i 

! 


I 

I 

I 

I 

I 

I 

I 

I 


Suppose  the  point  Q at  which  the  stresses  of  Figure  3 are 
to  be  determined  is  located  by  the  spherical  coordinates  shown 
in  Figure  4.  The  ranges  of  the  angles  0^,  (Jt-j  are  not  the  con- 


figure 4:  Alternative  spherical  coordinates. 

ventional  ones,  but  are  chosen  so  that  in  the  normal  (n  b^) 

-0-0 

plane  the  two-dimensional  form  of  the  stress  singularity  is 
recovered . 

In  the  transformation  of  coordinates,  equations  (8)  are 
replaced  by 

n^  ■ r^  cos0^  coS(^^ 

bo  - r^  sine^  (26) 

t©  " '*1  cos0^  sln^^ 
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Then  it  is  found  that 


s = r-|  cos0^  sin4»^  + 0(rp 

r sine  ■ sin6^  + O(r^)  (27) 

r cose  * cose^  cos(|>^  + O(r^) 


Therefore, 


tane, 
tane  = — ^ — 

(28)* 

r = r^XK  cose^ 


where 


X = cos^)^ 

. . ** 

C29) 

/ tane,  ' 

K * sgnCcose^)  / 1 + ( — ^— -^) 

The  definition  of  X in  equations  (29)  differs  from  that  in 
[6,9]  because  of  an  error  in  determining  the  relationship  be- 
tween the  elliptical  coordinates  and  the  spherical  coordinates 

The  terms  omitted  In  this  and  subsequent  equations  are  all 
one  order  of  magnitude  smaller  than  the  last  term  written. 

**The  symbol  for  curvature  Is  also  used  here  to  make  compari- 
sons with  [6,9]  easier.  There  should  be  no  confusion  because 
curvature  Is  contained  nowhere  In  Appendix  I. 
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near  the  crack  front.  Equation  (29)  Is  the  correct  form  to 
use. 

Equations  (28)  and  (29)  allow  one  to  compute  the  terms 
required  In  equations  (1).  Thus,  one  has 


tan0  = sgn(tane^)  /ic^-1 
sin0  = sgn(tan0^)  ^ 


cose  . i 


Sin  I 


cos  I 


sgn(s1n0^ ) 


sin  = sgn(sin0^ ) 


30  ^ 2-k  /nr 

cos  ^ ^ 


2^^K  / K-1 


(30) 


JL  0 = 1 sgn(s1n9^)  rZELZ 

/Ir  ^ /7r^  ^ 2XCOS0.J 

JL  rnc  9 - J 1 / IC-H 

/2r  ^ /StJ"  2Xcos0^ 

The  stress  field  obtained  by  substituting  equations  (30) 
Into  equations  (1)  and  using  equations  (18),  (19),  and  (27)  is 
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r 


° - (<*-<+  2)A  - — — (3k^  + k + 2)B 


= 


° - ^ (3ic^  + tc  - 2)A ^ (<^  - K - 2)8 


'^b  = 


/2r, 


/2r^ 


aj  = 2v  [ — ^ 2kM  - — ~ 2k2b] 


V?F7 


/Ir^ 


(31) 


nb 


/irj 


^ (k^  - K - 2)8  + (k*  - tc  + 2)A 


/TrJ 


nt 


^ 2k^8 


/Trj' 


tJ.  = — ^ 2tc*A 
bt  ^ 


where 


A 


1 /~k±L~ 

2 I ic|  * * 2Xcos9^  ’ 


8 


sgnCsine^)  y — ^ — 
2 1 K 1 ’ ixcose-i 


(32) 


Except  for  the  correction  of  X (equation  (29))  and  a more  pre- 
cise treatment  of  signs,  equations  (31)  agree  with  those  of 
[6.9]. 


The  complicated  form  of  equations  (31)  makes  it  difficult 
to  find,  except  by  numerical  analysis,  that  the  minimum  strain 
energy  density  occurs  when  > 0.  The  form  of  equations  (20) 
makes  the  result  obvious. 


-22- 


The  mixed  mode  loading  of  an  infinite  elastic  body  con- 
taining a flat  elliptical  crack  is  accomplished  as  in  Figure 
5 by  inserting  the  crack  at  an  angle  to  the  direction  of  a 


a 


Figure  5:  Elliptical  crack  at  an  angle  to 

a field  of  uniform  tension  in  an 
infinite  material . 

field  of  uniform  tension.  The  coordinates  of  point  on  the 
ellipse  are  given  in  terms  of  the  parametric  angle,  a,  by 


X = a cosa,  y = b sina 


(33) 


The  stress  intensity  factors  at  point  [7]  are 


a/a~ 

ETkT 


a/a~k'k^  sinB  cosB 

G fl  * 


r sinot  sinm 

*-(k^+vk*^)E(k)-vk'^K(k) 


k'  cosa  cosoj  1 
(k^-v)E(k)+vk'=^K(k)-‘ 


(34) 


(1-v)  a/aTk'k^  sing  cosS  [- 


k'  cosa  sing) 

(k'"  + vk'^)E(k)-vk'''^K(k) 


+ 


sina  COSO)  t 
CF-v)E(k)+vk'=^K(k)J 


where  a^  = b/l-k^cos^a 


and  K(k)  and  E(k)  are  complete  elliptic  integrals  of  the  first 
and  second  kinds,  respectively. 

These  stress  intensity  factors  can  be  substituted  into 
equation  (22),  and  the  fracture  angle  0*  determined  for  vari- 
ous angles  of  loading,  6 and  w,  crack  diameter  ratios,  b/a, 
and  positions  on  the  crack  front,  a.  The  point  of  the  crack 
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★ 

front  with  the  largest  value  of  the  minimum  S will  control 
the  load  required  to  cause  failure,  The  value  of  this 

critical  applied  stress  is  obtained  from  equation  (25). 

Several  results  in  [9]  are  modified  by  these  calculations. 

Table  9-1*  is  no  longer  applicable  since  the  minimum  value  of 

* * 

energy  density  occurs  in  the  normal  plane  . Table  1,  con- 
taining corrected  values  of  Table  9-2,  shows  the  direction  of 
minimum  strain  energy  density  corresponding  to  tensile  loading 
for  various  points  along  the  crack  front.  Also  given  are  the 
normalized  minimum  values  of  the  strain  energy  density  and  its 
distortional  (S^)  and  dilatational  (S^)  parts.  Corrected  val- 
ues for  Table  9-3  are  given  in  Table  2.  Since  the  correct 

it 

values  of  the  minima,  S , are  larger  than  those  reported  in 
[9],  the  failure  loads  will  be  smaller  than  in  [9J. 

Figure  6 shows  the  direction  of  crack  growth  at  various 
points  of  the  crack  front  (specified  by  a)  for  b/a  = 0.5.  The 
loads  are  applied  symmetrically  with  respect  to  the  xz-plane. 
As  a result,  the  ends  of  the  minor  axis  (a  = 90°)  are  under- 
going Mode  I and  III  deformation  only,  and  the  crack  grows 
straight  ahead  under  tension.  Other  values  of  the  fracture 
angle  are  given  as  functions  of  the  angle  of  inclination  of 
n 

The  labels.  Table  9-1,  Figure  9-1,  etc.,  will  be  used  to 
designate  Table  1,  Figure  1,  etc.,  in  reference  Labels 

without  the  9-prefix  refer  to  tables  and  figures  in  this  Re- 
port. 

**In  the  notation  of  [9],  » 0. 


25- 


Table  1.  Fracture  angles,  9*,  and  S*  for  b/a  = 0.5, 

V = 0.33,  o)=0,  and  3 = 60®.  (Tensile  load 
i ng) 


a 

* 

®o 

166S*/a^b 

16GS*/o2b 

16GS*/o2b 

0 

-47.55® 

0.4613 

0.1208 

0.3405 

15® 

-44.85® 

0.5015 

0.1535 

0.3480 

30® 

-38.50® 

0.5991 

0.2291 

0.3701 

45® 

-30.58® 

0.7144 

0.3144 

0.3999 

o 

o 

-21.56® 

0.8171 

0.3890 

0.4280 

75® 

-11.29® 

0.8870 

0.4394 

0.4476 

o 

o 

0 

0.9118 

0.4572 

0.4545 

Table  2.  Fracture  angles,  9*,  and  S*  for  b/a  * 0.5, 


- 0 

V = 0.33,  u=0,  and  3 
loading) 

= 60®. 
16GSj/o^ 

(Compressive 

a 

* 

16GS*/a2b 

b 16GS*/o^b 

0 

124.78® 

0.0884 

0.0778 

0.01056 

15® 

128.24® 

0.1168 

0.1095 

0.00730 

O 

O 

CO 

136.47® 

0.1851 

0.1823 

0.00272 

45® 

146.56® 

0.2648 

0.2643 

0.00059 

o 

o 

157.35® 

0.3355 

0.3355 

0.00006 

75® 

168.55® 

0.3835 

0.3835 

0.00000 

90“ 

180.00® 

0.4004 

0.4004 

0 

-26 


Fi^^CrTU AN6,LE^  />3  FUNCTlOf^^  C>F  OF  LO/^D'^i^ 

y/)icious  Fositio/v^  /}/vx>  h/a-o.s  F^OZl^:y-o. 


the  load  to  the  plane  of  the  crack.  As  expected,  the  frac- 
ture angle  1s  zero  when  the  load  Is  perpendicular  to  the  crack 
(&  = 90°).  Figure  6a  should  be  contrasted  to  Figure  9-4,  and 
6b  to  9-14.  The  curves  in  [9]  show  very  different  behavior 
near  the  ends  of  the  minor  axis.  Similar  statements  can  be 
made  about  Figure  7 In  comparison  to  Figures  9-5  and  9-15. 

It  shows  fracture  angles  for  a longer  elliptical  crack,  b/a 

= 0.1. 

Figure  8 shows  the  influence  of  the  length  of  the  ellipse 
on  fracture  angle  at  the  location,  a = 30°.  The  longer  el- 
lipses approach  tunnel  cracks  ttwo-dimenslonal ) which  for  o)=0 
are  under  Mode  I and  III  loading  only.  As  would  be  expected 
for  this  case,  the  fracture  angles  approach  zero  for  tensile 
loading.  Again,  the  curves  in  Figure  8 are  distinctly  dif- 
ferent from  Figures  9-6  and  9-7  for  the  smaller  values  of  b/a. 

Fracture  angles  for  several  values  of  the  angle,  w,  between 
the  xz-plane  and  the  plane  of  symmetry  of  the  loading  are  shown 
in  Figure  9.  These  angles  give  the  direction  in  which  fracture 
occurs  If  It  occurs  at  the  ends  of  the  minor  axis  as  Is  the 
case  for  Figures  13  and  14.  Figures  9-7  and  9-17  show  the  same 
curves  as  Figure  9 because  [9]  chose  to  put  equal  to  zero, 
which,  to  be  consistent  with  other  calculations  In  [9],  should 
have  been  zero  only  for  w=0  and  90°.  As  verified  In  this  Re- 
port, the  fracture  always  occurs  In  the  normal  plane,  and  there 
fore.  Inadvertently,  Figures  9-7  and  9-17  are  correct. 
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Figures  10,  11,  and  12  give  the  values  of  minimum  strain 
energy  density  for  various  values  of  the  parameters.  These 
curves  may  be  interpreted  as  giving 

S - F (35) 

* 

where  S is  the  minimum  strain  energy  density,  and  F,  as  shown 
in  the  figures,  is  a function  of  a,  3,  u,  v,  and  b/a. 

The  case  of  a narrow  ellipse  (b/a  = 0.1}  with  load  applied 
so  that  u)=0  is  shown  in  Figure  10.  The  curves  in  Figure  9-8 
and  9-18  are  qualitatively  similar,  but,  except  for  a=0  and 
90°,  there  are  large  differences.  The  case  of  tensile  loading 
shows  peaks  for  6 « 55°.  This  indicates  that  failure  is  more 
likely  to  occur  when  the  crack  and  load  make  this  angle  than 
when  they  are  perpendicular.  This  type  of  result,  already 
pointed  out  in  [8,9],  requires  rethinking  of  conventional  de- 
sign procedures  which  consider  the  worst  case  to  be  the  one 

* 

for  8 = 90°.  Compressive  loading  produces  peak  values  of  S 
at  8 2 45°  and  zero  values  at  8=0  and  90°.  This  implies  that 
extremely  large  loads  (which  would  cause  general  yielding  and/ 
or  crushing)  would  be  required  to  cause  failure  in  the  latter 
two  cases. 

* 

Figure  11  shows  the  effect  of  b/a  on  S . For  nearly  circu- 
lar cracks,  there  is  no  peak  under  tension  until  8 ” 90°. 
Therefore,  penny-shaped  cracks  are  in  the  most  detrimental  po- 
sition when  normal  to  the  load.  When  the  loads  are  compressive, 
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the  peaks  are  exactly  at  6 = 45®,  and  the  curves  are  symmetric 
about  6 = 45®  because  at  a = 90®,  only  Mode  I and  III  loading 
are  present  when  (jj=0.  Figures  9-9  and  9-19  are  Identical  to 
Figure  11  because  of  this  symmetry. 

As  discussed  In  connection  with  Figure  9,  Figures  9-10  and 

9-20  are  the  same  as  Figure  12  because  of  an  inconsistency  in 

[9]  which  Inadvertently  led  to  the  correct  result.  Figure  12 

* 

shows  values  of  S at  a = 90®  for  a long  ellipse  (b/a  = 0.1) 

and  several  orientations  of  the  load.  The  curve  for  uj=0  Is 

associated  with  Mode  I and  III  only,  while  to  = 90®  produces 

Mode  I and  II  only  at  a = 90®.  For  Intermediate  values  of  to, 

there  Is  a mixture  of  all  three  modes  of  crack  deformation. 

* 

The  curves  show  a decrease  In  S as  the  load  changes  from  a 
mixture  of  Modes  I and  III  to  one  of  Modes  I and  II. 

be  obtained  from  Figures  10,  11,  and  12 
C35)  and  (25).  Solving  for  the  failure 

(36) 
or 


where  F Is  the  function  of  a,  B,  u»,  v,  and  b/a  shown  In  Fig- 
ures 10,  11,  and  12.  The  form  used  in  equation  (37)  was  chosen 
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for  plotting  Figures  13  and  14. 

Figure  13  is  obtained  from  Figure  11  with  one  exception. 

For  most  cases,  it  was  found  that  the  largest  value  of  S 
occurs  at  a = 90“.  However,  for  nearly  penny-shaped  cracks, 
it  sometimes  occurs  at  a=0  instead.  Therefore,  the  case  of 
tensile  loading  shows  two  curves  for  b/a  = 1.0.  The  upper, 
solid  curve  shows  the  load  required  to  cause  failure  at 
a = 90®,  and  the  lower,  dashed  curve  gives  the  load  at  which 
failure  actually  occurs  (at  a=0).  It  is  only  in  the  inter- 
mediate range  of  3 that  the  dashed  curve  is  lower,  and  so  it 
is  not  shown  outside  this  range.  In  all  other  cases  plotted 
in  Figure  13,  the  failure  occurs  at  a = 90®  at  the  load  shown. 
It  is  difficult  to  compare  Figures  9-12  and  9-22  with  Figure 
13,  but,  by  symmetry,  they  should  give  the  same  result.  Note 
that  the  ordinates  in  £93  are  2/tt  times  those  in  Figures  13 
and  14. 

Failure  loads  for  the  various  mixtures  of  the  three  modes 
of  crack  deformation  are  shown  in  Figure  14.  This  was  obtained 
from  Figure  12  by  using  equation  (37).  In  all  cases,  the  fail- 
ure occurs  at  a = 90®  for  this  long  elliptical  crack  (b/a 
* 0.1).  It  can  be  verified  from  Figure  10  that  for  o)«0  the 
largest  value  of  S occurs  at  a = 90®. 

i 


-37- 


REFERENCES 


M.  L.  Williams,  "On  the  stress  distribution  at  the  base 
of  a stationary  crack".  Journal  of  Applied  Mechanics. 

Vol.  24.  pp.  109-114,  1957. 

2.  G.  C.  Sih,  "Stress  distribution  near  internal  crack  tips 
for  longitudinal  shear  problems".  Journal  of  Applied  Me- 
chanics. Vol.  32,  pp.  51-58,  1965. 

3.  R.  J.  Hartranft  and  G.  C.  Sih,  "The  use  of  eigenfunction 
expansions  in  the  general  solution  of  three-dimensional 
crack  problems".  Journal  of  Mathematics  and  Mechanics, 
Vol.  19,  pp.  123-138,  1969. 

4.  G.  R.  Irwin,  "Analysis  of  stresses  and  strains  near  the 
end  of  a crack  traversing  a plate".  Journal  of  Applied 
Mechanics . Vol.  24,  pp.  361-364,  1957. 

5.  I.  N.  Sneddon,  "The  distribution  of  stress  in  the  neigh- 
borhood of  a crack  in  an  elastic  solid".  Proceedings  of 
the  Royal  Society  of  London.  Series  A,  Vol.  187,  pp.  229- 
260,  1946. 

6.  M.  K.  Kassir  and  G.  C.  Sih,  Three-Dimensional  Crack  Prob- 
1 ems . Noordhoff  International  Publishing,  Leyden,  1975. 

7.  M.  K.  Kassir  and  G.  C.  Sih,  "Three-dimensional  stress 
distribution  around  an  elliptical  crack  under  arbitrary 
loadings".  Journal  of  Applied  Mechanics.  Vol.  33,  pp. 


601-611,  1966. 


-40- 


8.  G.  C.  Sih,  "A  special  theory  of  crack  propagation". 

Methods  of  Analysis  and  Solutions  of  Crack  Problems,  pp . 
XXI-XLV,  G.  C.  Sih,  editor,  Noordhoff  International  Pub- 
lishing, Leyden,  1973. 

9.  G.  C.  Sih  and  B.  C.  K,  Cha,  "A  fracture  criterion  for 
three-dimensional  crack  problems".  Engineering  Fracture 
Mechanics.  Vol.  6,  pp.  699-723,  1974. 

10.  G.  A.  Korn  and  T.  M.  Korn,  Mathematical  Handbook  for 
Scientists  and  Engineers,  pp.  561-584,  McGraw-Hill,  1968. 

11.  M.  E.  Kipp  and  G.  C.  Sih,  "The  strain  energy  density 
failure  criterion  applied  to  notched  elastic  solids". 
International  Journal  of  Solids  and  Structures.  Vol.  11, 
pp.  153-173,  1975. 

12.  R.  J.  Hartranft  and  G.  C.  Sih,  "Three-dimensional  growth 
characteristics  of  a plane  crack  subjected  to  concentrated 
forces".  Journal  of  Applied  Mechanics.  Vol.  41,  pp.  808- 
809,  1974. 


-41- 


DISTRIBUTION  LIST 


PART  I - GOVERNMENT 


Administrative  and  Liaison  Activities 

Chief  of  Naval  Research 
Department  of  the  Navy 
Arlington,  Virginia  22217 
Attn;  Code  439 
461 
444 

Director 

ONR  Branch  Office 
495  Sutimer  Street 
Boston,  Massachusetts  02210 

Director 

ONR  Branch  Office 
536  South  Clark  Street 
Chicago,  Illinois  60605 

Director 

Naval  Research  Laboratory 
Attn:  Library,  Code  2029  (ONRL) 
Washington,  O.C.  20390 


Headquarters 
U.S.  Marine  Corps 
Washington,  D.C.  20380 
Attn:  MC-AX-2  (2) 

MC-AX-4E2 
(3)  MC-AAP-3 

Conmander 

Naval  Air  Systems  Comnand 
Department  of  the  Navy 
Washington,  D.C.  20360 
Attn:  NAVAIR-320N 

NAVAIR-5302 
NAVAIR-531 
NAVAIR-5311J 
NAVAIR-5314 
NAVAIR-5325A 

Conmander 

Naval  Air  Development  Center 
Warminster,  Pennsylvania  18974 
Attn:  Mr.  M.  Schulman 
(6)  Crew  Systems  Department 


U.S.  Naval  Research  Laboratory 

Attn:  Technical  Information  Division 

Washington,  D.C.  20390  (6) 

Conmanding  Officer 
ONR  Branch  Office 
207  West  24th  Street 
New  York,  New  York  10011 

Director 

ONR  Branch  Office 
1030  E.  Green  Street 
Pasadena,  California  91101 

Defense  Documentation  Center 
Cameron  Station 
Alexandria,  Virginia  22314 

Navy 

Chief  of  Naval  Operations 
Department  of  the  Navy 
Washington,  O.C.  20305 
Attn:  N0P-05F 

N0P-506N 
N0P-098T 


Naval  Air  Development  Center 
Philadelphia,  Pennsylvania  19112 
Attn:  Dr.  E.  Handler 

Crew  Systems  Department 

Director 

Naval  Research  Laboratory 
Washington,  D.C.  20390 
Attn:  Code  8440  (Dr.  F.  Rosenthal) 

Commanding  Officer 
Naval  Ship  Research  and  Develop- 
ment Center 

Bethesda,  Maryland  20034 
Attn:  Code  745  (Mr.  B.  Whang) 

(12)  Conmanding  Officer 

Naval  Civil  Engineering  Labora- 
tory 

Port  Hueneme,  California  93041 
Attn:  Dr.  Warren  Shaw,  Structures 
Department 


Commander 

Naval  Safety  Center 
Naval  Air  Station 
Norfolk,  Virginia  23511 

Capt.  Channing  L.  Ewing  Mclusn 
Naval  Aerospace  Medical  Research 
Laboratory 

Michoud  Assembly  Facility 
New  Orleans,  Louisiana  70129 

Army 

Commanding  Officer 
U.S.  ^rmy  Aviation  Material 
Laboratories 

Fort  Eustis,  Virginia  23604 
Attn:  VDLEU-SS 

(Mr.  6.  T.  Singley,  III) 

Director 

U.S.  Army  Board  for  Aviation 
Accident  Research 
Fort  Rucker,  Alabama  36360 
Attn:  BAAR-PP  (Mr.  J.  Haley) 

Comnanding  Officer 

U.S.  Army  Research  Office  Durham 

Attn:  Mr.  J.  J.  Murray 

CRO-AA-IP 

Box  CM,  Duke  Station 
Durham,  North  Carolina  27706 

Air  Force 

Air  Force  Office  of  Scientific 
Research 

1400  Wilson  Boulevard 
Arlington,  Virginia  22209 
Attn:  Mechs.  Oiv. 

Commanding  Officer 
6571st  AMRL 

Holloman  Air  Force  Base, 

New  Mexico 

Attn:  Mr.  C.  C.  Gragg 
NASA 


Nr.  C.  Kubokawa 
N239-3 

NASA  Ames  Research  Center 
Moffett  Field,  California  94035 


Department  of  Transportation 

Mr.  R.  F.  Chandler 
Code  AC  119 

CAMI,  FAA  Aeronautical  Center 
P.O.  Box  25082 

Oklahoma  City,  Oklahoma  73123 
Mr.  H.  Spicer 

Federal  Aviation  Administration 
(DS-41) 

800  Independence  Avenue,  S.W. 
Washington,  D.C.  20590 

Mr.  H.  Daiutolo 

Federal  Aviation  Administration 
NAFEC  (NA541) 

Atlantic  City,  New  Jersey  08405 
Mr.  D.  Beyer 

Federal  Aviation  Administration 
(RD-730) 

800  Independence  Avenue,  S.W. 
Washington,  D.C.  20590 

Headquarters 
U.S.  Coast  Guard 
400  7th  Street,  S.W. 

Washington,  D.C.  20590 
Attn:  DAT/ 62 
EAE/63 
ENE-5/64 
IGS-1/61 
NMT-t/82 
OSR-2/73 

Mr.  J.  G.  Viner 

Protective  Systems  Group  (RS-12) 
Federal  Highway  Administration 
400  7th  Street,  S.W. 

Washington,  D.C.  20590 

Mr.  W.  H.  Collins 
Computer  Technology  Group  (DV-11) 
Federal  Highway  Administration 
400  7th  Street,  S.W. 

Washington,  D.C.  20590 

Mr.  Kenneth  Batcheller,  Chief 
Safety  Programs  Division  (RS-20) 
Federal  Railroad  Administration 
400  7th  Street,  S.W. 

Washington,  D.C.  20590 


Mr.  E.  Ward,  Chief 
Engineering  Research  and  Development 
Division  (RT-20) 

Federal  Railroad  Administration 
400  7th  Street,  S.W. 

Washington,  D.C.  20590 

Dr.  J.  A.  Edwards 

Associate  Administrator  for  Research 
and  Development 

National  Highway  Traffic  Safety 
Administration 
400  7th  Street,  S.W. 

Washington,  D.C.  20590 

Mr.  C.  D.  Ferguson 
Office  of  Crashworthiness  (41-40) 
National  Highway  Traffic  Safety 
Administration 
400  7th  Street,  S.W. 

Washington,  D.C.  20590 

Mr.  L.  L.  Bradford 

Office  of  Vehicle  Structures  Research 
(43-50) 

National  Highway  Traffic  Safety 
Administration 
400  7th  Street,  S.W. 

Washington,  D.C.  20590 

National  Transportation  Safety  Board 
800  Independence  Avenue,  S.W. 
Washington,  D.C.  20590 
Attn:  Mr.  A.  L.  Schmieg,  NS-10 

Mr.  H.  L.  Morgan,  NS- 20 

Mr.  B.  C.  Doyle,  NA-87 

Dr.  H.  E.  vonGierke 

Aerospace  Medical  Research 
Laboratory 

Aerospace  Medical  Division 
Air  Force  Systems  Command 
Wright-Patterson  Air  Force  Base, 

Ohio  45433 

Sal  Davis 

Fairchild  Industries,  Inc. 

Fairchild  Republic  Division 
Farmingdale,  New  York  11735 


Part  II  - CONTRACTORS  AND  OTHER 
TECHNICAL  COLLABORATORS 

Mr.  S.  P.  Desjardins 
Dynamic  Science 
1800  West  Deer  Valley  Drive 
Phoenix,  Arizona  85027 

Dr.  A.  A.  Ezra,  Chairman 
Department  of  Mechanical  Sciences 
and  Environmental  Engineering 
University  of  Denver 
Denver,  Colorado  80210 

Dr.  Albert  I.  King 
Bioengineering  Center 
Wayne  State  University 
Detroit,  Michigan  48202 

Dr.  R.  C.  DeHart,  Director 
Department  of  Structural  Research 
Southwest  Research  Institute 
P.O.  Drawer  28510 
San  Antonio,  Texas  78284 

Professor  J.  B.  Martin 
Division  of  Engineering 
Brown  University 
Providence,  Rhode  Island  02912 

Dr.  0.  L.  Tocher 

Boeing  Computer  Services  (72-80) 

P.O.  Box  24346 

Seattle,  Washington  98124 

Mr.  D.  G.  Harding,  Manager 
Survivability  Staff 
Boeing  Company  - Vertol  Division 
Philadelphia,  Pennsylvania  19142 

Dr.  H.  E.  Lindberg,  Manager 
Engineering  Mechanics  Program 
Stanford  Research  Institute 
Menlo  Park,  California  94025 

Mr.  John  W.  Freyler 
Beta  Industries,  Inc. 

2763  Culver  Avenue 
Dayton,  Ohio  45429 

Dr.  L.  E.  Hulbert,  Chief 
Advanced  Solid  Mechanics  Division 
505  King  Avenue 
Columbus,  Ohio  43201 


I 

I 


I 

I 

! 


Professor  George  Sih 
Department  of  Mechanics 
Lehigh  University 
Bethlehem,  Pennsylvania  18015 

Dr.  Harold  Liebowitz,  Dean 
School  of  Engineering  and 
Applied  Science 
George  Washington  University 
725  23rd  Street 
Washington,  D.C.  20006 


Dr.  Francis  Cozzarelli 
Division  of  Interdisciplinary  Studies 
and  Research 
School  of  Engineering 
State  University  of  New  York 
Buffalo,  New  York  14214 

Dr.  George  Herrmann 
Stanford  University 
Department  of  Applied  Mechanics 
Stanford,  California  94305 


I 

I 

I 

I 

I 


Professor  S.  B.  Dong 
University  of  California 
Department  of  Mechanics 
Los  Angeles,  California  90024 

Professor  A.  J.  Durelli 
Mechanics  Division 
The  Catholic  University  of 
America 

Washington,  D.C.  20017 

Professor  H.  H.  Bleich 
Department  of  Civil  Engineering 
Columbia  University 
Amsterdam  & 120th  Street 
New  York,  New  York  10027 


Professor  J.  D.  Achenbach 
Technological  Institute 
Northwestern  University 
Evanston,  Illinois  60201 

Professor  J.  Kempner 
Department  of  Aero.  Engrg.  and 
Applied  Mech. 

Polytechnic  Institute  of  Brooklyn 
333  Jay  Street 
Brooklyn,  New  York  11201 

Dr.  Nicholas  J.  Hoff 
Dept,  of  Aero,  and  Astro. 

Stanford  University 
Stanford,  California  94305 


I 

I 

I 

I 

I 

I 

I 

I 


Professor  A.  M.  Freudenthal 
George  Washington  University 
School  of  Engineering  and 
Applied  Science 
Washington,  D.C.  20006 

Professor  P.  G.  Hodge 
Department  of  Mechanics 
Illinois  Institute  of  Technology 
Chicago,  Illinois  60616 

Professor  D,  C.  Drucker 
Dean  of  Engineering 
University  of  Illinois 
Urbana,  Illinois  61801 

Professor  N.  M.  Newmark 
Department  of  Civil  Engineering 
University  of  Illinois 
Urbana,  Illinois  61801 

Library  (Code  0384) 

U.S.  Naval  Postgraduate  School 
Monterey,  California  93940 


Professor  Norman  Jones 
Massachusetts  Institute  of  Tech- 
nology 

Department  of  Naval  Architecture 
and  Marine  Engineering 
Cambridge,  Massachusetts  D2139 

Professor  Werner  Goldsmith 
Department  of  Mechanical 
Engineering 

Division  of  Applied  Mechanics 
University  of  California 
Berkeley,  California  94720 

Professor  W.  D.  Pi 1 key 
Department  of  Aerospace  Engineering 
University  of  Virginia 
Charlottesville,  Virginia  22903 

Dr.  H.  N.  Abramson 
Southwest  Research  Institute 
8500  Culebra  Road 
San  Antonio,  Texas  78206 


Mr.  John  Scowcroft 
Automobile  Manufactures  Association 
330  New  Center  Building 
I Detroit,  Michigan  48202 

Dr.  R.  D.  Young 

i Texas  Transportation  Institute 

Texas  A & M University 
College  Station,  Texas  77840 

I Professor  J.  A.  Collins 

Mechanical  Engineering  Department 
Arizona  State  University 
I Tempe,  Arizona  85281 


I 

I 

I 

I 

\ 

I 

I 


r 


Unci  ass i f i ed 

security  classification  of  this  rage  (Whtt  Dmts  Bnfrmd) 


REPORT  DOCUMENTATION  PAGE 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


2.  GOVT  ACCESSION  NO 


3.  RECIPIENT’S  CATALOG  NUMBER 


4.  TITLE 


Stress  Singularity  for  a Crack  with  any 
Arbitrarily  Curved  Fronts  ^ ^ 


tyre  of  rerort  a reriod  covered 


Technical  Rep 


6 PCRFoRSm^ORG.  REPORT  NUMBER 


7 AuTHORftJ 

J R.  J.^al 

' ^ • i 


( / ''' 


lartranft  ! 
G.  C./  Sih  /' 


J,  CONTRACT  OR  fiRAuT  NUypgR/'«> 

^ N0l)l4-76-C-de^94/  , 


9.  performing  organization  name  and  address 


Lehigh  UniyI 
InstH-ute  of  Fracture  & Solid  Mechanics, 
Bethlehem,  Pennsylvania  18015 


10.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  a WORK  UNIT  NUMBERS 


N0014-76-C-0094 


n.  CONTROLLING  OPEICE  NAME  ANO  AOORESS 

Department  of  the  Navy 
Office  of  Naval  Research 
Arlington,  Virginia  22217 


’?i._aE?OBT.JBA.TE 

Apr 


»9. 


PAGES 


t4-  monitoring  agency  name  a AOORESSff/ from  Controff/ng  Off/e») 


-43- 


15.  security  class,  (of  thio  eoport; 


15«.  DECLASSIFICATION/OOWNCRADING 
SCHEDULE 


16  distribution  statement  r®/ 


Approved  for  public  release;  distribution  unlimited. 


17  distribution  statement  (of  tho  sbotroct  ontorod  in  Block  70,  If  difforont  from  Roport) 


18  Supplementary  NOTES 


19.  KEY  WORDS  (Continuo  on  rovorom  oldo  If  nocooodry  ond  Idontity  by  block  numbor) 


strain  energy  density 
crack  front  stress  field 
\ crack  path 

^ree-dimensional  analysis 


^R>^T  ^Con(/nu«  on  rovoroo  oldo  If  noeooomry  ond  idontity  by  block  numbor) 

TO  i I 


20.  ABS 

nvestigate  the  possibility  of  a point  of  a crack  front 
propagating  in  a direction  not  in  the  normal  plane  of  the  crack 
front,  the  three-dimensional  form  of  the  crack  front  stress  field 
is  obtained.  To  simplify  comparisons  of  the  states  of  stress  at 
various  points  lying  on  a spherical  surface  centered  at  a point  o 
the  crack  front,  a local  spherical  coordinate  system  is  used.  It 
is  found  that  the  crack  propagation  will  be  from  each  point  of  th 


DD  /Sr,,  1473 


EDITION  OF  I NOV  (t  It  OBSOLETE 

// 


//nn  Uncusslflen /f 

(-fU  (/  f I SECURITY  CLASSIFICATION  OF  this  PAGE  (Whwi  0»»a 


*v*  «v  -a.a.p 


security  classification  of  This  RAOEfRImi  BaIa  Bn(Ai*0 


crack  front  in  a direction  lying  in  the  normal  plane. 

The  results  are  used  in  conjunction  with  the  strain  energy 
density  fracture  criterion  for  the  problem  of  an  elliptical  crack. 
The  plane  of  the  flat  elliptical  crack  makes  an  arbitrary  angle 
with  the  field  of  uniform  applied  tensile  stress.  Crack  growth 
directions  for  various  positions  along  the  crack  front  are  deter- 
mined, and  loads  required  for  fracture  for  various  angles  are 
obtai  ned .. 


Unclassified 


tlCURlTV  CVAStiriCATION  OF  THIS  RAOt0»AA  Oa<a 


INSTRUCTIONS  FOR  PREPARATION  OF  REPORT  DOCUMENTATION  PAGE 


^SpONSlpl^ITY  The  controlling  DoD  office  will  be  responsible  f.,r  compleuon  of  the  Report  Documentation  Page,  DD  Form  147.1.  m 
all  technical  reports  prepared  by  or  for  DoD  organisations. 

Ci*^^SlFJC A TiOT^.  Since  this  Report  Documentation  Page.  DI)  Form  1473,  is  used  in  preparing  announcements,  bibliographies,  and  data 
banks  it  should  be  unclassified  if  possible.  If  a classification  is  -.  quired,  identify  the  classified  items  on  the  page  by  the  soProDriate 
Svfrknnl  • * ^ 


COMPLETION  GUIDE 


General  Make  Blocks  1,  -1.  5.  6,  7.  11,  13,  IS,  and  16  agrcp  with  the  corresponding  information  on  the  report  cover  Leave 
Blocks  2 and  3 blank 

Block  1.  Report  Number.  Enter  the  unique  alphanumeric  report  number  shown  on  the  cover. 

Block  2.  Government  Accession  No.  Leave  Blank.  This  space  is  for  use  by  the  Defense  Oocumentstion  Center. 

Block  3,  Recipient's  Catalog  Number  Leave  blank.  This  space  is  for  the  use  of  the  report  recipient  to  assist  in  future 
retrieval  of  the  document. 

Blgtk  4.  Title  and  Subtitle.  Enter  the  title  in  all  capital  letters  exactly  as  it  appears  on  the  publication.  Titles  should  be 
unclassified  whenever  possible.  Write  out  the  En^ish  equivalent  for  Greek  letters  and  mathematical  symbols  in  the  title  (see 
“Abstracting  Scientific  end  Technical  Reports  of  Defense-sponsored  RDT  E.  “AD-667  000).  If  the  report  has  a subtitle,  this  subtitle 
should  follow  the  mam  title,  be  separated  by  a comma  or  semicolon  if  appropriate,  and  be  initially  capitalized  If  a publication  has  a 
title  in  a foreign  language,  tranalate  the  title  into  English  and  follow  the  English  translation  with  the  title  in  the  original  language 
Make  every  effort  to  simplify  the  title  before  publication  * 

Block  S,  Type  of  Report  and  Period  Covered.  Indicate  here  whether  report  is  inlenm,  final,  etc.,  and,  if  applicable,  in.  lusue 
dates  of  period  covered,  such  as  the  life  of  a contract  covered  in  a final  contraclor  report. 

Bl,ock  6.  Performing  Organization  Report  Number  Only  numbers  other  than  the  official  report  number  shown  in  Block  1 such 
as  series  numbers  for  in-house  reports  or  a contractor,  grantee  number  assigned  by  him.  will  be  placed  in  this  space  If  no  such  numbers 
are  used,  leave  this  space  blank. 

B'pck  7.  Authorfs).  Include  corresponding  information  from  the  report  covet.  Give  the  nameCsl  of  the  authoKsi  m conventional 
order  (for  example,  John  R.  Doe  or,  il  author  prefers,  J.  Robert  Doe,.  In  addition,  list  the  affiliation  of  an  author  if  it  differs  from  that 
of  the  performing  organization. 

Blocit  8.  Contract  or  Grant  Number's).  For  a contractor  or  grantee  report,  enter  the  complete  contract  or  grant  number/ .vy  under 
which  the  work  reported  was  accomplished.  Leave  blank  m in-house  reports. 

Block  9.  Performing  Organization  Name  and  Address.  For  in-house  reports  enter  the  name  and  address,  including  office  symbol, 
of  the  performing  activity.  For  contractor  or  grantee  reports  enter  the  name  and  address  ol  the  contractor  or  grantee  who  prepared  the 
report  and  identify  the  appropriate  corporate  division,  school,  laboratory,  etc  , of  the  author,  f.ist  city,  state,  and  ZIP  Code. 

Block  10.  Program  Element.  Project.  Task  Area,  and  Work  Unit  Numbers  Enter  here  the  number  code  from  the  applicable 
Department  of  Defense  form.  Such  as  the  DD  Form  1498,  "Research  and  Technology  Work  Unit  Summary”  or  the  DD  Form  1634 
"Research  and  Development  Planning  Summary."  which  identtfies  the  program  element,  project,  task  area,  and  work  unit  or  equivalent 
under  which  the  work  was  authorieed. 

PLock  U.  Controlling  Office  Name  and  Address.  Enter  the  full,  official  name  and  address,  including  office  symbol  of  the 
controlling  office.  (Equstes  fo  far^ding'sponfiorinti  agency.  For  de/inifion  see  OoD  Drreclive  S200.20.  **Distribvnon  Statements  on 
Technicaf  Documents.*') 

Block  12.  Report  Date.  Enter  here  the  day.  month,  and  year  or  month  and  year  as  shown  on  the  cover. 

Block  13.  Number  of  Pages.  Enter  the  total  number  of  pages. 

Block  14.  Monitoring  Agency  Name  and  Address  (it  differenf  from  Contro//ing  O/fue;.  For  use  when  the  controlling  or  funding 
office  does  not  directly  administer  a project,  contract,  or  grant,  but  delegates  the  administrative  responsibility  (o  another  organization. 

Blocks  15  & 15a.  Security  Classification  of  the  Report:  Declassificatlon.^Downgrading  Schedule  of  the  Report.  Enter  in  15 
the  hichest  classification  of  the  report.  If  appropriate,  enter  in  15a  the  declasaification/downgrading  achedule  of  the  report,  using  the 
abbreviations  for  declasaification/downgrading  schedules  listed  in  paragraph  4-207  of  DoD  5200. 1>R. 

Block  16»  Distribution  Statement  of  the  Report.  Insert  here  the  applicable  distribution  statement  of  the  report  from  DoD 
Directive*”S200.20,  ''Distribution  Statements  on  Technical  Documents." 

Block  17.  Distribution  Statement  (of  the  abstrat  ( tmfered  tn  Block  20.  if  different  from  the  distribution  statement  of  the  report). 
Insert  here  th^^plicable  distribution  statement  of  the  abstract  from  DoD  Directive  5200.20,  "Distribution  Statements  on  Technical  Doc- 
uments. " 


Block  lS.  Supplementary  Notes.  Enter  information  not  included  elsewhere  but  useful,  such  as:  Prepared  in  cooperation  with 
. . . TranaiaUon  of  for  by)  . . . Presented  at  conference  of  . . To  be  published  in  . . 

Block  19.  Key  Words.  Select  terms  or  short  phrases  that  identify  (he  principal  subjects  covered  in  the  report,  and  are 
Sufficiently  specific  and  precise  to  be  used  as  index  entries  for  cataloging,  conforming  to  standard  terminology.  The  DoD  "Thesaurus 
of  Engineering  and  Scientific  Tenna"  (TEST).  AD-fi72  000,  can  be  helpful. 

Block  20.  Abstract.  7>»e  abstract  should  be  a brief  fnof  fo  exi  ved  200  words)  factual  summary  of  the  most  significant  informa- 
tion contcTine^in  the  report.  If  possible,  the  abstract  of  a classified  report  should  be  unclassified  and  the  abatract  to  an  unclassified 
report  should  consist  of  publicly-  releasable  information.  If  the  report  contains  a significant  bibliography  or  literature  survey,  mention 
it  hare.  For  information  on  preparing  abstracts  see  "Abstracting  Scientific  and  Technical  Reports  of  Defense-Sponsored  RDT%E." 
AD-667  000. 


6S0  sso-ess 


